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ABSTRACT 

G : 

Context. Studies of the mass function (MF) of open clusters of different ages allow us to probe the efficiency with which brown 
dwarfs evaporate from clusters to populate the field. Surveys of older clusters (age > 100 Myr) are not affected so severely by several 

£ — ■ problems encountered in young clusters, such as intra-cluster extinction and large uncertainties in brown dwarf models. 

Aims. We present the results of a deep photometric survey to study the MF of the central region of the old open cluster Praesepe (age 

i i , 590^pjj Myr, distance 190+* g pc), down to the substellar regime. 

fV^ Methods. We performed an optical (riz and K-band) photometric survey of Praesepe using the Large Binocular Telescope Camera 

' covering an area of 0.59deg 2 in the cluster centre from i ~ 19.0 mag (~100M Jup ) down to a 5<x detection limit at i ~25.6 mag 

(~40Mj up ). The survey is approximately 95% complete at i = 23.8 mag and z = 22.0 mag (~55 Mj up ). 

i-C ' Results. We identify 59 cluster member candidates, of which 37 are substellar, by comparing with the predictions of a dusty atmo- 

sphere model. The MF of those candidates rises from the substellar boundary until ~67 Mj up and then declines. This is quite different 
JL 1 from the form inferred for other open clusters older than 50 Myr, but seems to be similar to those found in very young open clusters, 

the MFs of which peak at ~10Mj llp . Either Praesepe really does have a different MF from other clusters or they had similar initial 
MFs but a different dynamical evolution. Since most of the candidates are faint, we lack astrometric or spectroscopic follow-ups to 
test their memberships. However, the contaminations by field dwarfs, galaxies, or giants are found to have little effect on the shape of 

i i ■ MF and therefore the MF of 'real' cluster members should have similar characteristics. 



Key words, stars: brown dwarfs - stars: low-mass - stars: luminosity function, mass function - Galaxy: open clusters and associa- 
tions: individual: Praesepe 
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^^O , 1. Introduction ulations as a whole evolve, e.g., the efficiency with which BDs 

and low-mass stars evaporate from clust ers. These investig ations 

^ ■ The mass functions (MFs) of stellar and substellar populations have been carried out for the Hvad es dBouvier et al.l|2008l and 

have been determined from optical and near-infrared survey s for re f e rences therein) and for Praesepe (iBoudreault et al.ll2010l and 

, several open clusters at different ages, such as the Orion Nebula references therein) 

— _i ' Cluster, cr Orionis, p Ophiuchi, Taurus, IC 348, IC 2391, M35, „, „ , ■ , , ■ , • 

-.I. ™ ■ j ... TT j rr, ,, r , , , , The Praesepe open cluster has been surveyed extensively in 

• • , the Pleiades, and the Hyades. These MFs show clear heterogene- , , „ _ , , n-F , , . , , , J 

»> •» t tt irni u- u u n ju i . the past (cf. Tablellh, but only a few surveys have reached masses 

. _j ■ ity (see Fig. 110b, which may be partially caused by cluster evo- , , , , ,- ■ , , , , . s „ i t-.t^ 

J^j i Vj n — below the substellar limit (and then only just). Several BD can- 

r> . ' didates were detected in those surveys, some of which will be 

Studies of relatively old open clusters (age > 100 Myr) are re . examined in the t work . The su bstellar MF of Praesepe 

L « ' important for two particular reasons: first, they allow us to „„„,•„ ,, „. • 

. r, > r remains uncertain, 

study the intrinsic evolution of basic properties of brown dwarfs i 1 ■ , 

(BDs), e.g., luminosity and effective temperature, and to com- I Boudreault et al. | (|2010|) observed a significant difference be- 

pare the evolution with structural and atmospheric models; sec- tween the MFs of Praesepe and Hyades. While they found 
ond, we may investigate how the BD and low-mass star pop- H y ades MF has a maximum at -0.6 M dBouvier et al. | 

2008), the MF of Praesepe continues to rise from 0.8 M down 

to O.IMq. This is surprising, as both clusters share similar 

Send offprint requests to: W. Wang, e-mail: wangw@nao . cas . cn physical properties (ages, mass, metaUicity, and tidal radii). 

* The LBT is an international collaboration among institutions in Disagreement between the Praesepe and Hyades MFs could arise 

Germany, Italy and the United States. LBT Corporation partners £ . t . - a i , i • v i un £ j-xr 

T T , J ' „ .... ,, , „ v .. from variations in the clusters initial MFs, or from differences 

are: LBT Beteiligungsgesellschalt, Germany, representing the Max- . , . , . , , . r=r : -11 ' „ i ... , ... 

Planck Society, the Astrophysical Institute Potsdam, and Heidelberg ln thelr dynamical evolution dBastian et al. || 2010|) . Although dif- 

University; Istituto Nazionale di Astrofisica, Italy; The University of ferent bma ry fractions could cause the observed (system) MFs 

Arizona on behalf of the Arizona university system; The Ohio State to differ, there is no clear evidence of any variations in the bi- 

University, and The Research Corporation, on behalf of The University nary fractions from m easurements published in the literature 

of Notre Dame, University of Minnesota and University of Virginia. dBoudreault et al.ll2010h . 
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In this paper, we present a survey of the very low-mass 
star and substellar populations of Praesepe using the blue and 
red Large Binocular Cameras, extending down to hitherto unex- 
plored mass regimes (~40Mj up ). The main aims of our study are 
to search for new BDs and determine the MF of the Praesepe for 
a large coverage of the substellar regime. 

The candidate selection procedure, and mass and tempera- 
ture determination meth ods employed in t his study are similar 
to those adopted in Boudreault et aT] (120101) . However, we probe 
a lower mass regime and use an evolutionary model based on 
a dusty atmosphere instead of a combination of dust-free and 
dusty models. 



2. Observations and analysis 

2.1. Observations 

The observations presented in this paper were carried out with 
the Large Binocular Teles cope (LBT) located on Mount Graham, 
Arizona dHill etal.ll2 006). usin g the Large B nocular Cameras 
(LBCs, see lSpeziali etalJl2008h . The LBCs are two wide-field, 
high-throughput imaging cameras, namely Blue (LBCB) and 
Red (LBCR), located at the prime focus stations of the LBT. 
Each LBC has a wide field of view (23'x23'), with four CCD 
detectors of 2048x4608 pixels each, providing images with a 
sampling of 0.23"/pixel. 

The optical design and detectors of the two cameras are op- 
timized for different wavelength ranges: one for ultraviolet-blue 
wavelengths (320-500 nm, including the Bessel U, B, V and 
Sloan g and r bands), and one for the red-infrared bands (500- 
1000 nm, including the Sloan i, z and Fan Y bands). In the full 
binocular configuration, both cameras are available simultane- 
ously, and both point in the same direction of the sky, thus dou- 
bling the net efficiency of the LBT. 

To accomplish the entire survey of the inner region of 
Praesepe, we carried out three observing runs, in March 2008, 
December 2008, and February 2009. Table |2] summarizes the 
observations and Fig. 1 shows the areas surveyed. The total area 
covered is 0.59 deg 2 , about 1 percent of the cluster region. The 
transmission curves of the filters used in this survey is presented 
in Fig. 2, along with a synthetic spectrum of a brown dwarf with 

= 2300 K, log g = 4.5 [CGS], and solar metallicity (NextGen 
model). 



2.2. Reduction and astrometry 

The standard data reduction steps for the LBT data were per- 
formed using the IDL astronomy package and IRAF. The bias 
subtraction was executed on a nightly basis and for each CCD 
chip. To correct for pixel-to-pixel variations and global illumi- 
nation, master flat frames were created for the nights using twi- 
light exposures. For nights when no appropriate sky flat expo- 
sures were available, we used a master sky flat in the adjacent 
night. The individual images of a given field were registered 
and median combined, resulting in a combined science frame 
for each CCD, field and filter. To detect faint sources, we sub- 
tracted the strong background introduced by very bright stars3- 
We then used the IRAF task daofind to detect sources in the 
"clean" frames. The sources were extracted from the original sci- 
ence frames and instrumental magnitudes calculated using both 




Fig. 1. Spatial distribution of the LBT pointings 
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1 The area fraction affected by bright stars is less than 3% for most 
of CCD images, and is about 6% in the worst case. 
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Fig. 2. Transmission curves of the filters used in our survey 
compared to a synthetic spectrum of a brown dwarf with T e g 
= 2300 K, log g = 4.5, and solar metallicity (NextGen model). 
The transmission curves include the quantum efficiency of the 
detectors. 



aperture and point-spread function photometry with the IRAF 
tasks phot and allstar respectively. 

At this stage, weak fringes were stil l visible for several Y- 
band images. The method described by iBailer- Jones & M undt 
(2001) for removing fringes does not apply well in the present 
case, since the images are seriously affected by bright stars and 
no clean fringe images could be created. However, as the T-band 
images are ~1.5 dex shallower than expected, we decided not to 
use them in this study (although we still quote some statistics of 
the photometry below). 

An astrometric solution was achieved using the Sloan Digital 
Sky Survey (SDSS York et al. 2000) catalogue as a reference. 
The root mean square (rms) accuracy of our astrometric solution 
is 0.10-0.15 arcsec. As with other reduction procedures, astrom- 
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Table 1. Summary of previous photometric surveys in the Praesepe around the hydrogen-buming limit. The completeness limits 
correspond to a 10cr detection while those of the present work and Boudreault et al. (2010) correspond to 5<x detections. 



Authors 


Telescope / instrument 


Area 


BD 


Completeness limits 






(deg 2 ) 


candidates 


(mag) 


Hambly et al. (1995) 


COSMOS / POS & UKSTU 


19 





R F -20, I N -19 


Pinfield et al. (1997) 


TNT/WFC 


1.0 


-10° 


fl=21.5, /=20.0, Z=21.5 


Magazzu et al. (1998) 


INT / WFC 


0.22 


1 


fl=22.2, /=21.2 


Chappelle et al. (2005) 


INT / WFC 


2.6 


4 


7 C =21.3,Z=20.5 


Gonzalez-Garcfa et al. (2006) 


3.5 m CAHA / LAIC A & 5 m Hale / LFI 


0.33 


1 


/'=23.8, z'=23.3 


Boudreault et al. (2010) 


3.5 m CAHA / n2k & 2.2m La Silla / WFI 


3.1 


6 


/ c =23.2, 7=19.9, Ks=lS.6 


This work'' 


LBT / LBC 


0.59 


37 


r=24.1, /=25.6, z=24.7, 7=20.3 



° From Fig. 3 in Pinfield et al. (1997), about 10 of their 26 Praesepe member candidates have masses below 72 M lup . 
b Our work was complemented by the the JK S 3.5 m CAHA / £22k data from Boudreault et al. (2010). 



Table 2. Journal of the LBT/LBC observations 



Field 


RA 


DEC 


Date 


Filter 


(min) 


Seeing 
(arcsec) 


m itT 
(mag) 


PraeA 


08:40:53.76 


+ 19:52:41.9 


2009-02-28 


r 


114 


1.0 


25.8 








2009-02-28 


i 


42 


1.0 


25.8 








2009-02-28 


z 


54 


1.0 


24.8 








2009-02-28 


Y 


18 


1.0 


21.3 


PraeB 


08:39:14.23 


+ 19:52:41.9 


2009-02-28 


f 


54 


1.0 


25.1 








2008-03-06 


i 


42 


2.2 


25.7 








2009-02-28 


z 


54 


1.0 


25.2 








2008-03-04 


Y 


48 


1.4 


22.2 


PraeC 


08:39:14.36 


+ 19:27:18.0 


2008-03-07 


r 


90 


1.9 


25.6 








2008-03-07 


i 


30 


2.0 


25.8 








2008-03-07 


z 


30 


1.0 


25.1 








2009-02-28 


Y 


18 


1.0 


21.2 


PraeD 


08:40:53.63 


+ 19:27:18.0 


2008-12-29 


r 


99 


1.0 


25.8 








2008-12-29 


i 


42 


1.0 


25.8 








2008-12-29 


7 


39 


1.0 


24.8 








2008-12-29 


Y 


18 


1.0 


21.9 



" We also performed a shallower, 84-min pointing in r on 2008-12-30. 



etry was also performed separately for each CCD, to ensure that 
the solutions were as robust as possible. 

2.3. Photometric calibration 

To correct for Earth atmospheric absorption of the photometry, 
we calibrated the inferred data using the r, i, and z band values of 
Sloan Digital Sky Survey (SDSS) objects that were observed in 
the science fields. Zero point offsets were determined from the 
difference between the SDSS magnitudes and our instrumental 
magnitudes. Since these were obtained with objects in the same 
field of view for each science frame, we did not perform a colour 
or airmass correction when reducing our riz photometry. The er- 
ror introduced by this approximation is less than 0.05 mag. 

To calibrate our Y band photometry, we used our LBT 
; and z photometry and the Y band photometry from the 
United Kingdom Infrared Telescope Infrared Deep Sky Survey 
(UKIDSS, Lawrence et al. 2007). We found that the differences 
between Y band LBT magnitudes and Y band UKIDSS magni- 



tudes have a linear dependence on the ; - z colours, which can 
be described by the equation 

y UKIDSS ~ Y LB T,rmv = «() + « 1 X (/ - z), (1) 

On the basis of about 800 common objects between UKIDSS 
and our measurements, the ao,i coefficients were determined and 
the instrumental Y magnitudes were then transferred into the 
UKIDSS Y photometry system. For the same reasons as for our 
riz data, we did not (need to) perform a colour or airmass cor- 
rection for our 7-band photometry. 

The 5<x detection limits of our survey are ~ 25.6 mag and 
24.7 mag for the i and z bands, respectively. However, we do 
not expect all targets brighter than these limits to be able to be 
detected. We estimate the survey completeness by comparing 
the number of objects detected to the number predicted assum- 
ing a uniform three-dimensional spatial distribution of stars. As 
shown in Fig. [3] the number of detected sources in each band 
deviates from a log-normal relationship at bright and faint lim- 
its. From this, we estimate the completeness to our 5cr detec- 
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Fig. 3. Estimation of the completeness limit for our survey using the iz bands. The solid lines are the best linear fits before the turn 
offs. The vertical dashed lines are the 5<x detection limit and the vertical dotted lines are the saturation magnitudes. 



tion limit as 67.0% at i = 25.6 and 70.3% at z = 24.7 respec- 
tively, which corresponds to ~40Mj up assuming a dusty atmo- 
sphere. A similar estimation yields a completeness of 82.0% at 
r = 24.1 mag and 82.6% at Y — 20.3 mag, respectively. Our 
survey is approximately 95% complete at i = 23.8 mag and 
z = 22.0 mag. 

The relatively low completeness of our survey is possibly 
caused by the saturation of bright stars. For stars lying near to 
saturated stars, the photometric uncertainties are relatively large. 
A significant fraction of detected stars is then excluded because 
of their large photometric errors, which lowers the completeness. 
The total area seriously affected by bright stars in i band is 3-6 
percent. 

2.4. Candidate selection procedure 

The candidate selection introduced by Boudrea ult et al.1 ([2010) 
was adopted in the present work. Candidates were first selected 
based on the colour-magnitude diagram (CMD) using iz bands 
from our LBC observations^ A second selection was performed 



2 Our rY bands observations do not reach a similar stellar mass, hence 
are not used here. 



using a colour-colour diagram. While our rY bands observa- 
tions are not deep enough for our pre sent investigations, data 
from the near-IR photometric survey by iBoudreault et al.l (|20 10) 
- which fully covers our survey area, with a 5<x detection limit 
at ~55 Mj up in J and K s bands - was used instead for the second 
selection. In the third and final selection, we used the known 
distance to Praesepe to reject objects based on the discrepancy 
between their observed magnitude in J and the magnitude pre- 
dicted from the isochrones and our estimation of T e g. To be con- 
sidered as a cluster member, an object had to satisfy all three of 

these criteria. 

We use the evolutionary tracks of Chabri er et al. 
and the atmosphere models from lAllard et al 
- assuming a dusty atmosphere (the AMES-Dusty 
- to compute an isochrone for Praesepe using 
of 



2000) 



(2001) 



model) 
an age 

a+6.0 



59Q -i20 M y r dFossati et al.1 120081) . a distance of 



190^* gpc (Ivan Leeuwenl l2009f). and a sola r metallicity 



([Fe/H] = 0.038±0.039, |Friel & B oesgaardl Il992h. We ne glect 
the reddening [E(B -V) = 0.027+0.004 mag, |Taylor|[ 2006ll . The 
transmission curves we used for the filters for these calculations 
are plotted in Fig. [2] The effective temperature varies from 
500 K to 3900 K in steps of 100 K, while the gravity ranges 
from 4.0 dex to 6.0 dex in steps of 0.5 dex. 
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Fig. 4. CMD with i and z bands used in the first selection pro- 
cedure. As solid lines we show the isochrone computed from 
an evolutionary model with a dusty atmosphere (AMES-Dusty). 
The numbers indicate the masses (in Mj up ) on the model se- 
quence for various z magnitudes. The dashed lines delimit our 
selection band. 



2.4.1. Colour-magnitude diagram 

Candidates were first selected from our CMD by keeping all 
objects that are no more than 0.28 mag redder or bluer than 
the isochrones in all CMDs. This number accommodates errors 
in the magnitudes and uncertainties in the model isochrones. 
We also include the errors from the age estimate and distance 
to Praesepe. We additionally include objects brighter than the 
isochrones by 0.753 mag in order to include unresolved binaries. 
In Figure [4] we show the CMD where candidates were selected 
based on their values of z versus i-z. 

From a total of 44209 objects above the 5<x detection 
limit in i and z filters and below our saturation limit, 709 ob- 
jects are retained as candidate cluster members (98.4% are re- 
jected). Of these , 160 w ere detected in the NIR observations of 
Boud reault et al.l (HpiO). These objects are used in the selection 
process described below. 

2.4.2. Colour-colour diagram 

From this step, The candidate selection is based on both the op- 
tical iz data and the NIR JK S data. A candidate must be detected 
in every band. 

The second stage of the candidate selection involves retain- 
ing only those objects that lie within 0.28 mag of the isochrone 
in the colour-colour diagram. This value accommodates the 
photometric errors, uncertainties in the model isochrones, and 
the uncertainty in the age estimation of Praesepe. The colour- 
colour diagram with the selection limits is shown in Figure |5] 
where we also plot the theor etical colours of red gia nts (us- 
ing the atmosphere models of Hausc hTTdt et al.l[l999bl) and the 
theoretical colours of six galaxies with redshifts from to 2 
(Meisenheimer et al. 201 1). Neither the red giants nor the galax- 




1 2 3 4 5 6 7 
z - K s [mag] 



Fig. 5. Colour-colour diagram used in the second selection step. 
The solid line is the isochrone computed from an evolutionary 
model with a dusty atmosphere (the AMES-Dusty model; the 
masses in M for each z - K s colour are shifted for clarity). The 
dashed lines show our selection area. We also show the theo- 
retical colours of six galaxies as thin lines and the theoretical 
colours of red giants as thick lines. The six galaxies are two star- 
bursts, one Sab, one Sbc, and two ellipticals of 5.5 and 15 Gyr, 
with redshifts from z=0 to z-2 in steps of 0.25 (evolution not 
considered). We assume that all red giants have a mass of 5 M Q , 
0.5 < log g < 2.5 and 2000 K < r eff < 6000 K. 



ies are expected to be a significant source of contamination; most 
of the low redshift galaxies that were not automatically discarded 
during PSF photometry with full-width-half-maxima (FWHMs) 
broader than the average stellar FWHM by -30%, were rejected 
by means of the visual inspection of individual cluster mem- 
ber candidates after selection procedures. Of the 160 objects se- 
lected in the first step, 88 are kept here. 

2.4.3. Observed magnitude vs. predicted magnitude 

As indicated in Section [3~2l our determinations of T e s and mass 
are based on the spectral energy distribution of each object, so 
are independent of the assumed distance. The membership status 
of an object can therefore be assessed by comparing its observed 
magnitude in a band with its magnitude predicted from its T e g 
and Praesepe's isochrone (which assumes a distance). The pre- 
dicted magnitude of a background contaminant would be lower 
(brighter) than its observed magnitude and higher (fainter) for a 
foreground contaminant. To avoid removing unresolved binaries 
that are real members of the cluster, we keep all objects with a 
computed magnitude of up to 0.753 mag brighter than the ob- 
served magnitude. We also take into account photometric errors 
and uncertainties in the age and distance of Praesepe. This selec- 
tion procedure is illustrated in Figure |6] 

Of the 88 objects selected through CMDs and colour-colour 
diagrams in the first two steps, 74 are retained here. After this 
step, we perform a direct visual inspection of the images to re- 
ject resolved galaxies and spurious detections. This inspection 



6 



Wang, W. et al.: Mass function of Praesepe from deep LBT observations 



40 



50 



Mass (M Jup ) 
60 70 



80 



90 



a 
a 
o 
2 



Q 

w 
> 
a 
&a 
m 
CD 

o 



• • • 



1500 



8000 
r.„ (K) 



2500 



Fig. 6. Difference between the observed J magnitude and the 
model J magnitude computed from the derived mass and T e ff, 
as a function of T e ff. The vertical line marks the location of 
L0 dwarfs. The dotted line (at -0.753 mag) represents the offset 
due to the possible presence of unresolved binaries, the dashed- 
dotted lines represent the error in the magnitude determination, 
and the long-dashed lines represent the uncertainties in the age 
and distance of Praesepe. The horizontal solid line just traces 
zero. 



removes 15 objects from the photometric selection, of which 
seven are possibly galaxies and eight are false detections. 

The remaining objects constitute our final cluster member 
candidates, shown as large dots in Fig. [7] Those selected using 
only the iz photometry amount to 709, and are presented in Fig. [7] 
as small dots. We note that employing NIR JK S data helps us to 
remove a significant fraction of contaminations. At each mass 
bin, we calculate the number of stars removed as a result of in- 
cluding JK S data in the selection, and use this to estimate the 
number of stars in the final mass bin where JK S data are unavail- 
able. 

ISchmidt etail d2010l) and lWest etail d2011l) investigated the 
colors of L and M dwarfs, respectively, for every spectral types 
using the SDSS and Two Micron All Sky Survey (2MASS) cata- 
logues. The observed colour ranges are consistent with our com- 
puted colour ranges for M and early L dwarfs. However, the ob- 
se rved i - J color r ange fo r individual L dwarfs as shown in Fig. 2 
in ISchmidt et alj ((2010) is 3.5-6.0 mag, broader and slightly 
redder than our model ranges (the dashes lines in our Fig. 5). 
W e note that the Hamm er spectral-typing procedure employe d 
bv ISchmidt et al l d2010l) , was developed bv lCovev et"ail (120071) . 
who sought to optimize only for K and M dwarfs. For L dwarfs, 
the uncertainties should be ~2 subclasses or larger. In addi- 
tion, the Hammer is designed for solar-metalli city dwarfs, whil e 
Praesepe has a metallicity of 0.27+0. lOdex dPace et alJl2 008). 
hence will be prone to significant, systematic errors (ICovev et all 
|2007|) . We therefore keep the use of our model colours alone for 
candidate selection for a homogeneou s study. We found t hat if 
we use the observed color range from Sch midt et al.l (120101) in- 
stead, the observed rise in the mass function around 60 Mj up (cf. 
Section 4) remains - our main conclusion remains unchanged. 
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Fig. 7. Same as Fig. |4] but with large dots standing for the 59 
final candidates, and small dots for the 709 candidates that pass 
the first selection procedures. 



3. Results 

3.1. Selected photometric candidates 

We find that 59 photometric candidates survive the selection pro- 
cedures (based on isochrones assuming dusty atmospheres), a 
density of about 100 objects per square degree. Details of all 
photometric candidates are listed in Table 3. The identification 
number (ID) of a candidate is defined according to the field in 
which it was found and a sequential number for that field. The 
last column, J mo dei, is the predicted J magnitude based on pho- 
tometric determination of T e ff and mass. 

Our survey concentrates on the substellar regime. Saturation 
occurs at ~18 mag in z band, corresponding to ~100Mj up . 
Therefore, most of t he low mass candida t es discovered i n pre- 
vious surveys (e.g. iPinfield et ail 1 19971 lHamblv et ail 1 19951) 
are saturated in o ur LBT images. Only a few faint brown 
dwarfs clas sified bvlPinfield et al.l (1 19971) . iGonzalez-Garcfa et all 
(2006), and Boudreault et al. (2010J) are rediscovered in the cur- 
rent survey (cf. Table 4). These objects are M5-9 dwarfs, ac - 
cording to the photometric relations given bv lWest et al1 (f2008). 

Some of our targets were previously identified as cluster 
members but rejected by our selection procedures or visual in- 
spection. For example, eleven of the 150 Boudreault et al. (2010) 
candidates are detected in our LBT survey (the rest are mostly 
saturated). However, nine of them are rejected in this work (cf. 
Table 5). Among them, seven are rejected on the basis of the z 
vs. i—z CMD selection, because they are bluer than the isochrone 
area. Another one is obviously not a point-like source in the 
LBT image, and another is rejected because its observed J mag- 
nitude is inconsistent with its model-predicted magnitude. The 
remaining two targets are confirmed to be cluster dwarf stars. 
As our current work e mployed more photometric bands than 
Boudreau lt et al.l (1201 Ol) . it is unsurprising that our selection is 
more conservative. 

Most of our candidates are in the substellar regime, and other 
than the five targets listed in Table 4, no other accurate pho- 
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Table 4. Previously discovered low-mass Praesepe candidates 
found in our survey. 



Table 6. Number density of the cluster member candidates. 



ID RA(J2000) 



A16 
A17 
C12 

D12 
D08 



DEC(J2000) 

(° "') 



Alternative name 



41:24.916 
41:26.483 
39:47.778 



8:40:53.559 
8:40:39.292 



Ref." 



19:57:26.19 NGC2632 PHJ20 [1] 

19:51:59.73 J084126.5+ 195200 [2] 

19:28:03.10 Praesepe017 [3] 

NGC2632PHJ11 [1] 

19:41:00.05 PraesepeOOl [3] 

19:28:39.49 J084039.3+ 192840 [2] 



Re fe: [1] iPi nfield et alj (199 
m lBoudreault et al.U20iaT ~ 



[2] iGonzalez-Garcfa et alj 120061 : 



logM 


T ff 


r 

" model 


Nfcand.) 


N(cont.) 


Fraction 


(M Jup ) 


(K) 


(mag) 


(deg" 2 ) 


(deg" 2 ) 




1.625 


1412 


23.29 


13 


195 


1500% 


1.675 


1692 


20.28 


27 


41 


150% 


1.725 


1839 


19.27 


18 


3.3 


15% 


1.775 


1981 


18.87 


33 


4.8 


13% 


1.825 


2244 


18.01 


38 


2.3 


5.7% 


1.875 


2361 


17.71 


17 


3.2 


16% 


1.925 


2479 


17.43 


8.4 


1.7 


17% 


1.975 


2668 


17.01 


8.4 


1.6 


16% 



tometric observations are available from past epochs. This pre- 
cludes using proper motions as a means of cluster membership 
assessment at this time. 

3.2. Photometrically-derived masses and effective 
temperatures 

From the evolutionary tracks and atmosphere models described 
previously, we obtained the magnitudes m.A from the average flux 
of a star in a specific band A using the equation 

m A = -2.5 log/ A + c A , (2) 

where niA is the magnitude observed in a given passband, /a 
the flux received on Earth in this passband, and ca is a con- 
stant that remains to be determined. The flux /a was obtained 
using the total transmission function of the passband for a given 
filter, convolved with the quantum efficiency of the CCDs (we 
assumed that the telescope and instrumental throughput is flat 
over each passband). To transform our optical i and z band mag- 
nitudes from the models to the Johnson photometric system, we 
assumed all magnitudes niA to be equal to 0.03 when /a is the av- 
erage Vega flux received on Earth. The constant ca f or each pass- 
band i s then determined using the Vega flux from IColina et al.l 
(fl992h . The fitted values for these constants are i = -22.180 
and z - -22.706 mag. The values computed by Boudr eault et al.l 
(2010) for J and K s were -23.687 and -25.908 mag, respec- 
tively. 

The masses an d effective temperatures wer e estim ated in the 
way described by Bou dreault & Bailer- Jonesl (|2OO90. For the 
faintest objects where only i and z are available, the colour i-z is 
used to compute masses and T e g. There are different sources of 
errors for the estimation of the mass and T e s, including the pho- 
ton noise, the photometric calibration, the least squares fitting 
(imperfect model), and the uncertainties in the age of and dis- 
tance to Praesepe. The latter two are the most significant errors 
and give the uncertainties of 0.008 M G and 263 K for a 0.05 M 
substellar object ( r eff = 1690K ), O.O1OM and 260K for 
a 0.06 M Q substellar object ( r eff = 1990K ), and 0.008 M 
and 201 K for a 0.072 M object at the hydrogen burning limit 
(r eff =2293K). 



3 We first normalized the energy distribution of each object to the en- 
ergy distribution of the model using the J filter. The energy distribution 
was then fitted via a least squares fit of the model magnitudes to the 
measured ones. 



3.3. Contamination by non-members 

As mentioned above, the three main sources of contamination 
are background red giants, unresolved galaxies, and Galactic M 
and L dwarfs. Red giants contaminate the high mass end of this 
study, as seen in the i — J vs. z - K s diagram, hence can be 
ignored. Although some types of galaxies have similar colours 
to Praesepe cluster members more massive than 60 Mj up , these 
low-redshift galaxies are in general extended sources and there- 
fore easily rejected by our visual inspection. Among the 74 can- 
didates that passed our selection procedures, we identify four 
as galaxies on the basis of their LBT images. Other possible 
sources are field L d warfs and high redsh ift quasars (for in- 
stance at redshift z ~6: ICaballero e t al. 2008). However, because 
these quasars have spectral energy distributions similar to mid-T 
dwarfs, whereas our faintest candidates have colours of early L 
dwarfs, and given that they are rare (0.25 quasars at 5 .5 < z < 6.5 
in a 0.59 deg 2 survey, Stem et al. 2007), the MF should not be 
affected by quasar contamination. 

The contaminatio n by field dwarfs is not negligible. 
ICaballero et al.l (120081) identified possible field dwarf contami- 
nants covering spectral types from M3 to T8 from the literature 
and presented the spatial density in the solar neighbourhood in 
their Table 3. From this, the spatial density of field dwarfs in the 
vicinity of Praesepe can be easily inferred, given the Galactic 
latitude of Praesepe of b — +32.5 deg and its distance of 190pc, 
assuming an exponential decrease for stellar density perpendic- 
ular to the Galactic disk with a scale height of 500 pc. The ab- 
solute J band magnitude range constrained by our selection pro- 
cedures is ~ +2 mag, as shown in Fig. [6] corresponding to a 
certain distance interval and a survey 'volume', which is defined 
by the product of survey area and depth. We calculated the num- 
ber of contaminants by multiplying the survey volume by spatial 
density at each mass bin. The result is shown in Table [6] The 
first column gives the central value of log M in each interval, 
while the second and third columns present the corresponding 
T e g and / mo dei values at that mass. The fourth column is the num- 
ber density of cluster member candidates (also shown Fig. |8]as 
filled triangles) after applying all corrections (except for in the 
lowest two bins, where contaminations are too high to be cor- 
rected). The final two columns give the number density of field 
dwarf contaminants and corresponding fraction. We found that 
the field dwarf contaminants do not affect the MF shape. The 
contaminat ion is significant for 7 m d e i ^ 20 mag. At this mag- 
nitude, the Boudreault et al.l (120101) / band has a completeness 
of 88%. Therefore, below this magnitude, the mass functions we 
derived are probably upper limits. 
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We conclude that the various contaminants are either neg- 
ligible, or do not affect the MF shape in the range that we can 
investigate quantitatively, i.e., from about 53 to 94M Jup . 

4. Mass function of very low mass and substellar 
population of Praesepe 

The mass function, ^(logioM), we present here is the total num- 
ber of objects per square degree in each logarithmic mass inter- 
val logioM to logioM + 0.1. Since we do not make any correc- 
tions for binaries, we compute here a system MF. 

As our candidates have been selected only from their pho- 
tometric properties, cluster membership confirmation via spec- 
troscopy is desirable. However, these observations will not be 
feasible in the near future because of the faintness of our can- 
didates. The following discussion is therefore based on the as- 
sumption that the MF of candidates is similar to that of 'real' 
cluster members. The assumption is possibly valid because our 
derived MF is consistent with that given by Boudreault et al. 
(2010) in the common mass bin, and the contamination by field 
dwarfs, giants, and galaxies should not affect the shape of the 
MF significantly. 

To account for the survey detection efficiency, we use a sim- 
ple simulation. For example, to calculate the detection efficiency 
of candidate A01 in i band, we select a bright but unsaturated 
star in the i band CCD image in which A01 resides, scale it 
down to the magnitude of A01 (i.e., i = 20.29), and randomly 
cast this 'fake' star in the CCD image 100 times. We then search 
and re-measure the 'fake' star with our procedures. The detec- 
tion efficiency is the fraction of 'fake' stars that have been re- 
discovered. We run this test for each filter (izJK s ) and multiply 
the detection efficiencies together (as we need a detection in ev- 
ery filter) to evaluate the overall detection efficiency for each 
candidate. This detection efficiency (or recovery rate) is about 
90% for the brightest candidates and drops very quickly to 10% 
for the faintest candidates. 

We mentioned in Section [2] that our optical photometry 
reaches lower masses than the NIR photometry that we used. 
To compute the MF of Praesepe to the lowest mass bin reached 
without optical data, we first computed a MF0 using only the op- 
tical iz photometry. This 'MF' is presented in Fig.|8]as filled dots. 
We computed a second MF from the list of candidates that pass 
the three selection crite ria and are also detected in the survey of 
Boudreault et al. (2010) in the NIR J and K s bands (presented on 
Fig. [8] as filled triangles). For each mass bin, we computed the 
number of objects removed as a result of adding the J and K s fil- 
ters to our selection process and mass determination procedure 
(plotted as a function of mass in Fig. [8] top panel). We fitted 
a linear function to estimate the number of objects that would 
be removed ;/ we had additional J and K s photometry to 40- 
45 Mj up , which is our lowest mass data point in Fig. [8] However, 
as shown in Table 6, the contamination in the two lowest mass 
bins is so overwhelming that the MFs in these two bins can only 
be regarded as upper limits and are no longer discussed in the 
paper. 

Our derived MF (presented in Fig. 8, 9, & 10) shows a rise 
from 105 Mj up to 67 Mj up and then a turn-over at ~67 Mj. This 
turn-over occurs well above the 5<x of either iz bands or JK S 
bands (e.g. at ~67 Mj up , i ~22) and we note that we have cor- 
rected the incompleteness of our survey and field dwarf contam- 
inations. We therefore believe that this feature is genuine. This 
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4 Note that this is an inaccurate 'MF', because of serious contamina- 
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Fig. 8. Lower panel. Mass function based on (a) the LBT i and 
z photometry (filled dots) after detection efficiency corrections, 
and (b) the LBT photom etry plus the J and K s photometry from 
iBoudreault et al.1 d2010h (dotted open triangles: the original MF; 
filled triangles: the MF after detection efficiency corrections and 
removing contaminants). The selection and mass calibration are 
based on dusty atmospheres. Error bars are Poissonian arising 
from the number of objects in each bin, except for the lowest- 
mass bin, for which the error bar is from the linear fit in the top 
panel. The vertical thin dotted lines are the mass limits at which 
detector saturation occurs in the i and z-bands. The vertical thin 
long dashed line is the mass at the 5cr detection limit of our 
optical LBT data while the thick short dashe d line is the mass 
at the 5<x detection limit of the NIR data of IBoudreault et ail 
(2010). Top panel. Difference of the log of the number of objects 
(in each mass bin) between the MF computed using the optical 
iz da ta and the MF compute d using both iz and NIR JK S data 
from Boudreault et al . (2010). The dotted line is a linear fit. 



is the first time a clear rise in the substellar MF in an old open 
cluster has been observed. 

The MF of Praesepe near the hydrogen-burning limit 
was previously obtain ed in several studies. However, only 
Boudreau lt et al.l (1201 Oh provide a common mass range for com- 
parison, as shown in Fig. [9] In the first substellar mass bin 
(~80Mj up ) we see that, both surveys give consistent results 
within their error bars. However, fo r the second bin a t ~70 Mj UD , 
our MF is much higher than that of Boudrea ult et al.l (120101) ; the 
discrepancy is smaller when considering the MF from that work 
using the dusty atmosphere (open dots), which is still a reason- 
able model for such low mass stars (~M9/L0). This may indicate 
that s ome faint candidates are missing in the Boudreault et al.l 
(2010) survey, as these authors did not make any corrections for 
the detection efficiencies. 

We emphasize that our LBT survey cover s the very central 
0.59 d eg 2 of Praesepe, while the Q2k survey by Boudreault et al. 

covers a much wider area (~3.1 deg 2 ). If no significant 
candidates are missing in the Q2k survey, this discrepancy may 
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suggest that the very low mass cluster members are mostly con- 
centrated in the cluster centre, in contrast to what is expected 
from a 'dynamical evaporation' of brown dwarf in open clus- 
ters. The basic idea of dynamical evaporation is that lower mass 
stars in a cluster have higher speeds according to equipartition 
of energy, so are able to move higher in the gravitational po- 
tential well of the cluster. Hence the fraction of low mass stars 
should increase with increasing distance from t he cluster centre- 
By co mparing the Praesepe and Hyades MFs, iBoudreault et al.l 
(2010) concluded that Praesepe might have been less affected by 
dynamical evolution. 

Owing to its large distance and old age, no other 
published MF determination of Praesepe has reached 
masses below 70Mj up . We therefore compare our results 
with those from other clusters in Fig. [1 01 This includes 
IC 239 1 from IBoudreault & Bailer- JonesT (120091) . ONC 
from Hillenbrand & Carpenter (2000), cr Orionis from 
dCaballero et all 120071; iBihain et alT "l2009). and the Hyades 
from Bouvier et al J (I2008h . The MF of Praesepe is quite dif- 
ferent from both IC2391 (age of -50 Myr) and the Hyades 
(~625 Myr). Either the 'dynamical evaporation' does not have 
(or has not yet had) the same effect on these three clusters, or 
they had different initial mass functions. Another possibility 
is that Praesepe has a different binary fraction. Employing 
different cluster member selection criteria may also account for 
the observed MF discrepancies among clusters. Further studies 
are necessary to clarify these points. 

The continuing rise of the MF into the substellar regime 
that we observe has also been observed in young clusters (a s 
shown in Fig. [10), especially in cr Orionis (Bihainet al. 2009), 
Trapezium (turn-over at ~10-2 0Mt, ip . | Muenchet aTl 12002). 
pOph (MF rising to ~10M Jup , Mars h et al.ll2010l) . and in the 
very low luminosity young cluster in S 106, where the MF in- 
crease s or at least remains flat down to ~10Mj up dOasa et al. I 
2006). If we assumed a universal IMF, then it seems that the 
substellar MF of Praesepe has not evolved significantly since the 
cluster formed. 



5. Conclusions 

We have carried out the deepest survey to date of the old 
open cluster Praesepe, covering the central 0.59 deg 2 in the rizY 
bands. The survey probed a mass range from ~100 to 40Mj llp 
at 5cr detection limit, with which we have derived the very low 
mass and substellar mass function of this cluster. 

We compared our optical /z-bands data, c ombined with the 
Q2k N IR (J and K s ) band observations from Boudreault et al.l 
(2010), with theoretical loci of cluster members based on a dusty 
atmosphere (the AMES-Dusty model), to select cluster mem- 
ber candidates. Our final sample comprises 59 photometric can- 
didates. We estimate that the contamination by field dwarfs is 
about 15%, and that this does not affect the shape of MF. The 
contamination by galaxies and red giants is believed to be negli- 
gible. About two thirds of our cluster members have theoretical 
masses below the hydrogen-burning limit at 0.072 M , and are 
therefore brown dwarf candidates. We emphasize that to claim 
cluster memberships for the candidates, follow-up astrometric or 
spectroscopic observations are required. However, given that the 
candidates are generally faint and these observations are very 
time-consuming, none of them has yet been confirmed in this 
way. The discussion in this contribution therefore refers to the 
mass function of photometric cluster member candidates. 

The mass function we have inferred for the central 0.59 deg 2 
of Praesepe is consistent with that inferred for a wider area 
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Fig. 9. MF of Praesepe based on our survey LBT iz and 
Q2k JK S photometry (triangles), compared with that from 
Boudreau lt et al.1 (1201 Ol) (open dots assuming a dusty atmo- 
sphere and filled dots assuming a dust-free atmosphere). Error 
bars are Poissonian arising from the finite number of objects ob- 
served in each bin, except for the last bin, for which the error bar 
is derived from the linear fit. The vertical thin dotted line is the 
mass limit above which detector saturation occurs in the i and 
z-bands. The vertical thin long-dashed line is the mass at the 5<x 
detection limit of our optical LBT data. 



by IBoudreault et al.l (120101) at a mass just below the substel- 
lar boundary, but deviates by ~0.6dex in the next lowest mass 
bin, which may indicate that there is either a significant num- 
ber of objects missing in the Boudreault et al. 2010 survey, 
or a higher concentration of substellar objects in the centre of 
Praesepe (as the Boudreault et al. survey is at a larger cluster ra- 
dius). The latter possibility suggests that the dynamical evolution 
of very low mass stars is inefficient in this cluster, as proposed by 
Boudreau lt et al.l (|2010) for explaining the discrepancy between 
the Praesepe MF and Hyades MF. 

The steady rise in the Praesepe MF down to ~70Mj up and 
a turn-over there were unexpected for this old cluster. Such a 
significant peak has never been observed in any other cluster 
older than 50 Myr, but has been observed in several very young 
open clusters such as cr Orionis or clusters in star-forming re- 
gions (e.g., Trapezium). This suggests that the dynamical inter- 
actions in Praesepe have very little effect on MFs, if we assume 
there is a universal initial MF. 
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ID RA(J2000) DEC(J2000) i err z err J err K s err T eff M J moiA 

(mag) (mag) (mag) (mag) (K) (M Jup ) (mag) 

A01 8:40:12.599 19:56:50.33 20.29 0.03 19.01 0.06 18.14 0.05 17.01 0.05 3023 

A02 8:40:18.893 19:57:07.40 20.17 0.02 18.74 0.06 17.75 0.03 16.77 0.04 2729 

A03 8:40:20.624 19:43:40.72 23.42 0.06 21.75 0.10 19.55 0.11 17.88 0.10 1921 

A04 8:40:24.901 19:57:15.52 22.93 0.05 21.19 0.03 19.25 0.18 17.53 0.08 1940 

A05 8:40:25.044 19:41:44.56 19.23 0.02 18.40 0.02 17.37 0.02 16.47 0.03 3093 

A06 8:40:25.773 19:57:08.28 19.84 0.03 18.79 0.04 17.63 0.04 16.72 0.04 2821 

A07 8:40:25.981 19:57:07.02 19.72 0.03 18.69 0.04 17.55 0.04 16.63 0.04 2900 

A08 8:40:27.459 19:43:19.69 22.06 0.08 20.66 0.08 18.60 0.04 16.96 0.04 1997 

A09 8:40:32.282 19:43:53.72 25.33 0.15 23.23 0.07 20.28 0.20 18.84 0.24 1855 

A10 8:40:33.639 19:53:08.76 23.94 0.03 22.41 0.02 19.77 0.06 18.37 0.08 1902 

All 8:40:35.313 19:44:54.81 19.72 0.01 18.80 0.02 16.66 0.01 15.83 0.02 2351 

A12 8:40:44.975 20:00:24.80 23.43 0.02 21.83 0.02 19.09 0.02 17.56 0.04 1842 

A13 8:40:49.194 19:59:12.88 23.19 0.02 21.59 0.03 19.31 0.04 18.41 0.08 2040 

A14 8:41:08.239 19:47:08.93 23.13 0.09 21.55 0.04 19.05 0.02 17.69 0.03 1905 

A15 8:41:15.526 19:44:11.43 20.49 0.02 19.20 0.01 17.00 0.07 16.13 0.03 2169 

A16 8:41:24.916 19:57:26.19 19.11 0.01 18.12 0.06 16.13 0.01 15.25 0.01 2356 

A17 8:41:26.483 19:51:59.73 21.38 0.01 20.13 0.02 18.24 0.02 17.34 0.02 2294 

B01 8:38:39.176 19:42:54.85 20.21 0.01 18.99 0.01 17.30 0.03 16.26 0.02 2373 

B02 8:38:59.208 20:02:32.63 22.99 0.03 21.31 0.01 19.29 0.12 17.95 0.10 2050 

B03 8:39:06.931 19:51:00.19 19.28 0.01 18.26 0.09 16.98 0.01 16.13 0.02 2712 

B04 8:39:28.285 19:44:01.63 19.62 0.01 18.69 0.09 17.64 0.04 16.66 0.03 3093 

B05 8:39:34.775 19:43:06.65 22.61 0.05 21.28 0.03 19.11 0.07 17.77 0.09 2060 

B06 8:39:37.010 19:52:30.06 22.20 0.02 20.72 0.03 18.32 0.05 16.68 0.03 1892 

B07 8:39:38.172 20:01:17.52 19.32 0.01 18.49 0.01 17.04 0.02 16.06 0.02 2674 

B08 8:39:51.841 19:50:42.94 22.63 0.02 20.90 0.03 18.31 0.05 16.84 0.03 1867 

B09 8:39:56.429 19:47:31.41 22.96 0.03 21.17 0.04 19.36 0.09 17.95 0.11 2112 

C01 8:38:34.501 19:31:08.83 19.27 0.01 18.44 0.01 16.81 0.09 16.13 0.02 2590 

C02 8:38:48.106 19:31:43.17 23.21 0.06 21.28 0.02 19.73 0.10 18.14 0.12 2181 

C03 8:38:49.426 19:16:43.34 23.50 0.03 21.74 0.02 18.98 0.06 17.34 0.05 1808 

C04 8:38:59.835 19:31:52.00 21.71 0.01 20.37 0.01 18.62 0.04 17.59 0.07 2310 

C05 8:39:00.756 19:37:38.70 19.08 0.01 18.21 0.01 16.70 0.01 15.77 0.01 2626 

C06 8:39:02.279 19:37:38.47 22.32 0.01 20.85 0.01 18.72 0.07 17.88 0.08 2159 

C07 8:39:11.773 19:36:35.63 21.53 0.01 20.27 0.01 18.61 0.10 17.44 0.06 2381 

C08 8:39:22.309 19:38:37.41 22.06 0.01 20.35 0.01 18.42 0.10 16.96 0.04 2062 

C09 8:39:35.180 19:34:13.22 19.97 0.01 18.94 0.01 16.50 0.10 16.41 0.03 2160 

C10 8:39:35.196 19:33:50.49 20.21 0.01 19.16 0.01 17.51 0.01 16.47 0.02 2462 

Cll 8:39:36.737 19:34:55.25 24.59 0.12 22.83 0.03 20.16 0.18 18.43 0.13 1786 

C12 8:39:47.778 19:28:03.10 20.15 0.01 18.80 0.01 16.90 0.01 15.82 0.01 2217 

C13 8:39:54.504 19:20:09.00 23.06 0.03 21.55 0.01 19.16 0.07 18.16 0.11 2007 

D01 8:40:10.438 19:24:07.42 19.39 0.01 18.49 0.02 17.08 0.01 16.25 0.02 2683 

D02 8:40:13.215 19:27:00.51 24.64 0.07 22.92 0.04 20.14 0.19 18.06 0.11 1703 

D03 8:40:13.829 19:26:55.82 24.33 0.05 22.58 0.05 20.35 0.24 18.65 0.20 1896 

D04 8:40:20.994 19:38:46.21 24.11 0.14 22.15 0.07 19.72 0.17 18.10 0.12 1848 

D05 8:40:21.561 19:38:51.23 23.88 0.15 21.81 0.06 19.31 0.12 17.40 0.06 1760 

D06 8:40:25.698 19:36:25.68 19.52 0.01 18.56 0.04 16.89 0.01 16.10 0.02 2504 

D07 8:40:26.220 19:37:54.58 23.70 0.08 22.28 0.05 19.71 0.13 18.44 0.16 1947 

D08 8:40:39.292 19:28:39.49 22.04 0.02 20.44 0.01 18.07 0.02 16.76 0.02 1962 

D09 8:40:43.705 19:29:52.89 25.40 0.17 23.39 0.10 21.10 0.38 18.66 0.14 1705 

D10 8:40:50.371 19:29:46.29 22.55 0.05 21.22 0.04 18.88 0.05 17.20 0.04 1898 

Dll 8:40:53.452 19:19:13.98 23.02 0.03 21.33 0.03 18.82 0.04 17.10 0.03 1847 

D12 8:40:53.559 19:41:00.05 19.85 0.02 18.60 0.03 16.82 0.05 15.61 0.02 2320 

D13 8:40:56.008 19:25:31.85 24.62 0.10 22.55 0.06 20.01 0.12 17.67 0.06 1682 

D14 8:41:00.974 19:32:00.88 23.01 0.08 21.23 0.06 19.08 0.06 17.81 0.07 2015 

D15 8:41:02.632 19:22:13.66 23.66 0.04 22.09 0.04 20.22 0.16 19.35 0.25 2216 

D16 8:41:06.062 19:27:47.71 24.41 0.08 22.88 0.06 20.26 0.17 19.41 0.30 1903 

D18 8:41:15.768 19:28:15.81 19.79 0.01 18.83 0.02 17.35 0.01 16.49 0.02 2588 

D19 8:41:17.045 19:28:13.46 22.84 0.02 21.41 0.02 19.44 0.08 18.38 0.12 2184 

D20 8:41:32.924 19:32:13.40 24.01 0.03 22.40 0.04 19.52 0.09 17.86 0.07 1807 

D21 8:41:36.404 19:25:06.45 20.71 0.01 19.60 0.01 18.01 0.03 17.24 0.05 2513 
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Table 5. The 9 Praesepe cluster candidates identifie d by Boudreault et"aT] (2010h that were rejected in this work by colour or 
brightness. The IDs in the first column are taken from Boudreault et al J (120101)" The fourth object is not a point-like source in the 
LBT image, thus no optical photometry measurement is available. 



ID RA(J2000) DEC(J2000) i z h J K s M r cff i raodcl 

(mag) (mag) (mag) (mag) (mag) (M Q ) (K) (mag) 



005 


8:41:08.50 


+ 19:54:02.0 


20.19 


18.70 


19.02 


16.58 


15.39 


0.088 


2636 


17.06 


007 


8:39:39.56 


+ 19:47:54.3 


18.93 


17.80 


17.95 


16.10 


15.07 


0.104 


2860 


16.58 


009 


8:39:55.84 


+ 19:53:14.3 


18.05 


17.78 


20.29 


17.50 


16.54 


0.081 


2520 


17.32 


018 


8:39:42.79 


+ 19:35:48.2 






18.27 


16.20 


15.21 


0.097 


2782 


16.78 


022 


8:41:04.20 


+ 19:31:27.8 


20.12 


19.15 


18.89 


16.67 


15.75 


0.092 


2702 


16.97 


901 


8:39:59.45 


+19:43:37.4 


18.05 


18.30 


19.09 


17.16 


16.41 


0.084 


2574 


17.20 


902 


8:39:23.72 


+ 19:52:01.8 


20.41 


19.46 


20.15 


17.77 


16.88 


0.073 


2348 


17.72 


903 


8:40:00.20 


+ 19:30:27.0 


19.59 


19.01 


19.74 


17.50 


16.62 


0.076 


2409 


17.57 


914 


8:38:52.02 


+ 19:35:05.3 


19.48 


18.69 


19.12 


17.25 


16.35 


0.085 


2591 


17.16 



